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Abstract. Binary stars produce an array of dramatic astrophysical phenomena. They 
allow us to probe stellar structure, nuclear physics, and gravitational wave physics. 
They also produce the powerful supernovae that allow us to measure the scale of the 
universe. Despite their importance and ubiquity, many questions remain unanswered 
as to how the star formation process produces the wide array of stellar multiples that 
we observe. A complete model for binary formation encompasses three main compo- 
nents. We must know: (1) the primordial population of systems, (2) the influence of 
the dynamical processes that reshape this distribution as stars form and natal clusters 
disperse, and (3) the role of binary stellar evolution. In this article I review the most 
prominent theories for binary formation: turbulent core fragmentation, disk fragmenta- 
tion, and competitive accretion. I argue that turbulent core fragmentation at all masses, 
with disk fragmentation added in at the upper end of the mass spectrum, might explain 
the trend towards increasing multiplicity for higher mass stars. In addition, I provide 
a brief overview of the observational statistics and of some of the important processes 
that modify the primordial distribution of stellar orbits. 



1. Introduction 



The study of binarity dates back at least 250 years to Mitchell (1767 1, who pointed 
out the abundance of double stars. Many tout the 21st century as the era of precision 
cosmology, but we are also entering the era of precision in more classic astronomy - 
understanding the ordering of the stars. Great technological advances including the use 
of coronography ( [Hinkley et al.|2007] ) and adaptive optics ( [Lafreniere et al.|2009] l have 



brought stellar companions into better focus. 

Theories of star formation have been influenced by our heliocentric bias. Until 
recently, our models have reflected that we hve next to a relatively low mass star (ap- 



parently) without a stellar or substellar companion. The seminal work of Duquennoy 



& Mayor] ( jl991( ) on the multiplicity of sun-like stars showed that our sun is somewhat 



unusual in being alone. Indeed if we consider planets as companions, very few stars are 
truly isolated ( |Youdinj2011p . Recent work by Ra ghavan et al.| ( |20l0p has confirmed that 
roughly half of sun-like stars are in binary or higher order multiple systems. Although 
the statistics are incomplete for stars of other masses, it seems clear that there is a strong 
correlation between stellar mass and multiplicity (Lada 2006; Mason et al.|2009 1. The 



M stars, the most common stars in the galaxy, are usually in single systems. As primary 
mass increases to just a few solar masses, binaries become ubiquitous, even as the stars 
themselves become quite rare. The data are consistent with all O and B stars having 
at least one stellar companion at birth. The statistics on higher order multiplicity are 
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still incomplete. It remains difficult to detect relatively close, extreme mass ratio pairs, 
particularly for massive primaries. Upcoming searches with next generation telescopes 
such as ALMA and JWST may improve these statistics. 

Armed with this new data, it is time to reconsider our models of single and mul- 
tiple star formation. To be sure, conclusions drawn from the current data require many 
caveats. We do not have a complete picture of companions at all separations and mass 
ratios. Moreover, given the diverse nature of these systems, it is unlikely that one mech- 
anism can explain them all. This review aims to provide a broad, though necessarily 
incomplete, picture of the state of the field. 

I begin by reviewing theories for binary formation, and proceed to discuss the cur- 
rent observational constraints on both the primordial and field populations. Tightening 
our focus to systems that might serve as progenitors of compact binaries requires the 
consideration of another set of processes: those that alter the orbital configurations of 
the birth population. In this review I will touch only on mechanisms that operate during 
the main phase of star formation, leaving late stage evolution to other chapters in the 
proceedings. 

A detailed examination of the properties of observed binaries is beyond the scope 
of this article. For an update to date observational census of low and high mass multiple 
systems see Raghavan et al. (2010 1 (low mass) and Sana & Evans (2010) (high mass). 
For a discussion of higher order multiples, see also |Tokovinin '{ 2008 1. 



2. Theories of Binary and Multiple Star Formation 



There are several prominent theories of binary and multiple star formation in the lit- 
erature. The formation mechanisms can be divided roughly into three categories (1) 
core mediated processes, (2) disk mediated processes, and (3) few-body dynamical 
processes. In this review, I neglect channels of binary formation that occur after star 
formation has ended. There is strong observational evidence that most binary formation 
is concurrent with star formation: the binarity fraction for pre-main sequence stars is 
higher than that for main sequence stars (Mathieu|1994 Kraus et al.|2011 1. 



2.1. Core Mediated Theories 

The most common mechanisms in the literature fall into the first category, and are 
described alternately as prompt or early fragmentation, or more specifically as core 
or turbulent fragmentation poyle||1953[ [Tsuribe & Inutsuka|[T999at [Boss et al.||2000[ 



Padoan & Nordlund||2002| Fisher||2002 l. The earliest fragmentation scenario is that 
of |Hoy le (1953) who argued for opacity-limited hierarchical fragmentation, whereby 
an unstable cloud keeps fragmenting at the local Jeans mass, which decreases as the 
cloud collapses to higher densities, until the smallest fragments become optically thick 
and can no longer cool efficiently. Although this mechanism is inconsistent with the 
observed stellar mass function (for a thorough review see Fisher^2002) , it gave birth to 
more modern theories of core fragmentation. Another classic binary formation theory is 
known as core fission, where a core deforms through a series of equilibrium ellipsoidal 
figures. Cores contract on too short a timescale (compared to their viscous timescale) 
to evolve in this manner (see also |Tohlme| ( |2002] ). 

In the subsequent generation of "core" fragmentation scenarios, a collapsing, bound 
gas clump was thought to fragment into two or more objects depending on the ratio of 



The Formation of Binaries 



3 



thermal to gravitational energy, a - 5cIr/GM, and the ratio of rotational to gravita- 
tional energy, /? ^ ffR^/3GMm where Cs is the sound speed, R is the core radius, M 
the core mass, and Q. the orbital frequency of the core. Inherent in this parameterization 
is the assumption that any dynamically important magnetic fields have already diffused 
out (or can be included simply as an extra effective pressure in the calculation of a). 
When the product a/? < 0.12, cores were initially thought unstable (Inutsuka & Miyama] 
1992[ ), although this model has the somewhat counterintuitive result that slowly rotating 



(small y6) virialized cores remain prone to fragmentation. More recent work has shown 
this criterion to be incomplete, requiring also that a < 0.5 ( [Tsuribe & Inutsuka. 1999b^ ). 
Yet this model neglects many thermal effects as the gas collapses to higher densities. 
The source of angular momentum setting yS in the model above was thought clas- 



sically to be galactic shear (e.g. Bodenheimer 1995 and references therein), but the 



ubiquitous turbulence in the interstellar medium could also be responsible for the ob- 
served line-width size relations interpreted as evidence of core rotation ( |Fisher|[2004 
Goodman et al.|I993| l. 



This identification between rotation and turbulence has led to yet another class of 
models, appropriately termed turbulent fragmentation. The turbulence scenario sug- 
gests that non-linear perturbations expected in a turbulent cloud can cause a sub-region 
within a core to become over dense and collapse more rapidly than the free-fall timescale 
of the background core, thereby leading to the production of a secondary condensation 
within the bound core. Alternatively, turbulent motions can lead to the formation of 
filamentary structures, which then fragment into multiple objects. These stars are pre- 
sumed to accrete from their natal core mostly independently. Galactic shear may help 
stir the turbulence in molecular clouds, although stellar feedback can also be important 
dFall et al.|2010l ). 



Recent work by Offner et al. ( 2009 1 and Offner et al. ( 20 1 1 have demonstrated that 



the turbulent fragmentation mechanism is the dominant source of binaries in radiation 
hydrodynamic simulations of low mass (M < M©) star formation. These simulations 
also show rapid orbital migration over the first few tens of thousands of years follow- 
ing formation. Although the binaries were bom on scales of ~ O.lpc, their separations 
had shrunk to less than 1000 AU after ~ lO'^ years. The final separation distribution 
produced by turbulent fragmentation remains unknown: simulations which follow the 
interaction of the gas with stars and the radiation field remain too computationally ex- 
pensive to run for millions of years. However, one expects the mass ratio distribution 
to mirror the IMF if the cores accrete onto their own stars relatively unmolested by the 
binary partner. 

Note that there is another class of star formation models dubbed turbulent frag- 
mentation or "gravoturbulent fragmentation" which refers to the fragmentation of a 
molecular cloud on larger scales into many bound objects which form their own stellar 
systems (Ballesteros-Paredes et al. 2007] ). The turbulence sets the mass-scale of the 



cores, and their turbulent structures, but for the purposes of multiple formation we are 
concerned with the subsequent fragmentation of the turbulent clump, not the process 
that generated it. 

2.2. Disk Driven Formation 

Models for binary and multiple formation at a slightly later evolutionary state involve 
protostellar disks. In the disk fragmentation model (Adams et al.|1989; Laughlin & Bo- 



denheimer|1994 Bonnell|1994 Bonnell & Bate|1994b I stars form sufficiently massive 
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protostellar disks that become susceptible to gravitational instabilities. If disks become 
unstable, and the gas can cool efficiently, they can fragment to produce one or more 
coplanar companions that accrete from the parent disk, and depending on the formation 
epoch, the primary's natal cloud. 

Disk fragmentation has received much attention from numericists in recent years. 
Prior work on this topic has been severely limited by two assumptions. The first is that 
disks can be modeled in isolation. Disks modeled in the absence of ongoing accretion 
are much less likely to become unstable. Ongoing accretion at the outer edge of the 
disk is nearly always necessary to drive instability (Kratter et al.|2008 2010b, ,Kratter| 
&Murray-Clay|2011| l. 



The second related assumption was that disks were always low in mass relative 
to the central star. Certainly observed T-Tauri disks tend to be lower in mass, but as 
observations improve we are beginning to observe more massive disks at early epochs 
of formation (Andrews & Willia ms|2007 Andrews et al.(2009[ |. Together these two no- 
tions led most previous authors to conclude that disk fragmentation might occur rarely, 
and could only produce very unequal mass ratio objects, which was inconsistent with 
binary statistics. 

However, if disks are more massive and undergo continued accretion, simulations 
suggest that fragmentation may be more common, and that secondaries born in the disk 
can accrete quickly to catch up to the primary in mass ( Bonnell & Bate|1994a Kratter 



et al.|2 010a', though see Hanawa et al.|2010 1. Together with observations of more mas- 



sive disks, this suggests that the disk fragmentation scenario can in fact produce a wide 
range of mass ratios. We predict that disk fragmentation is most relevant for systems 
with primaries more massive than our sun, however recent observations of young clus- 



ters indicate disk fragmentation might be relevant at lower masses (Kraus et al.^2011 



In particular, Kratter et al. (20 10a I showed that two dimensionless parameters, which 
measure the infall rate compared to the sound speed = M/{c]/G), hereafter normal- 
ized accretion rate) and the infall rate compared to the orbital time (F = M/{Mtd^, 
hereafter normalized rotation period), can accurately predict when disk fragmentation 
occurs. We show a compilation of these results in Figure 1. 

Disk fragmentation is a promising mechanism for the formation of relatively close 
binaries, because the separations at birth are less than lOOOAU, and often within a 
hundred AU. Due to interactions with their natal circumstellar or circumbinary disks 
(see significan t migration can also occur, much as in the case of planets ( jGoldreich 



& Tremaine|T979||Artymowicz & Lubow|1996| l. 



Disk fragmentation has also been shown to produce higher order, hierarchical mul- 
tiples (Krumholz et al.|2007 2009 ; Stamatellos et al.|201 1 1. Such systems may be ideal 



progenitors to compact binaries (see Perets & Kratter 201 1, submitted). 

As with the turbulent fragmentation model, disk fragmentation has yet to provide 
strong predictions for the separation distribution of binaries (or multiples). In addition, 
the mass ratio distribution is unclear. While a range of outcomes are possible, the mass 
ratios will be regulated by the relative timescale of fragmentation with respect to the 
accretion time from the core. In systems where fragmentation happens early, stars are 
likely to approach equal mass as they share infalling material; in contrast, systems that 
only fragment at the tail end of the star formation process will likely remain at disparate 



masses. The latter case might lead to the production of brown-dwarf companions (Sta- 



matellos & Whitworth.2009^ . 



The Formation of Binaries 



5 




Figure 1. Illustration o f different disk outc omes as a function of dimensionless 
accretion rates defined in Kratter et al. ( 2010a i. The disks on the right hand side have 
infall properties characteristic of more massive stars (M > 2Mq), while those on the 
left hand side better resemble low mass stars. The increasing propensity towards disk 
fragmentation at higher masses might explain the increased multiplicity fraction for 
stars of types earher than A. 
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2.3. Multiple Formation in a dynamically evolving cluster 



A final mechanism for multiple formation is that put forth in a series of papers includ- 
ing patel ( pOOOl ), |Bate et aLl ( |2003] l, [Bonnell et al.|p004l ), and |aarketaL] ( pD08| ). These 
modelers suggest that star clusters form through a dynamic process dubbed competi- 
tive accretion. In this scenario, a turbulent molecular cloud forms many roughly Jeans 
mass clumps that interact with one another and compete for mass from the background 
cloud. In this environment there are multiple mechanisms for binary formation, in- 
cluding the disk methods mentioned above, as well as dynamical interactions between 
three or more bodies which we discuss below. In addition, dynamical friction between 
the protostars and background gas is thought to tighten wider binaries. What distin- 
guishes this scenario from those listed above is that a stellar system does not accrete 
from a fixed, marginally bound core, but rather migrates around in a cloud collecting 
previously unbound material. 

A model put forth by Goodwin & Kroupa ( 2005| ) relies on a combination of the 
core fragmentation and dynamical evolution hypothesis to reproduce the high multi- 
plicity fraction of young stars in clusters, compared to their field star counterparts. 
Although they do not specify the precise mechanism for fragmentation, these authors 
imagine that most cores either break up into 2 or 3-i- individual stars. Through dynami- 
cal interaction, the higher order multiples evolve into single star systems, tight binaries, 
and stable hierarchical systems. They posit that different cluster environments (which 
have different rates of close stellar encounters) will naturally evolve to show different 
multiplicity fractions. 

A related mechanism is disk capture. In the capture scenario, a massive disk in- 
creases the gravitational cross-section for capturing a nearby star into a bound orbit. 
Although initial investigations found this process to be highly inefficient around low 



mass stars at reasonable stellar densities ( Clarke & Pringle 1991) , more recent work 
suggests that the process might be better suited to massive stars with correspondingly 
more massive disks (Moeckel & Bal ly||2007 1. However, because this process remains 



inefficient even in dense clustered environments, it is unlikely to be the dominant for- 
mation channel. 

With these models in hand, we now turn to mechanisms that reshape the primordial 
distributions. 



3. From Birth to Death: Mechanisms for Evolution 

Two types of interaction are likely to modify the birth population of binaries and multi- 
ples while the star formation process is ongoing: (1) accretion disk interactions, and (2) 
dynamical interactions in few body systems. We address the role of each of these pro- 
cesses in making compact binary progenitors. As noted above, molecular gas may also 
play a role in shrinking the separations of newborn binaries, however this could well 
be considered part of the core fragmentation process rather than a separate mechanism 
since this process operates before the stars have well defined orbital parameters. 

3.1. Disk Interactions 

Binaries born via any mechanism, and especially disk fragmentation can reside within a 
circumbinary disk of material as they accrete from the background cloud. Gravitational 
torques can expel the circumbinary disk carrying away a large fraction of the orbital 
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angular momentum ( |Syer & Clarke] 1995 ; Artymowicz & Lubow| 1996 1. The extent 
to which this process shrinks binary orbits is uncertain, and depends on the mass ratio 
between the stars and the disk. Some simulations suggest that once the disk is expelled 
to roughly twice the orbital separation of the binary, orbital evolution stalls, leaving a 
stable circular binary ("Arty mowicz & Lubow|1996 ^. This mechanism can likely bring 
binaries down to separations of lO's of AU, and possible much closer. 

Sometimes, both stars in a binary retain their own circumstellar disks. In this 
case, tidal interactions between each of the two circumstellar disks allows for transfer 
of angular momentum from the disks to the orbit, expanding it as in the case of the 
Earth-Moon system where the rotation of the disk / Earth is faster than that of the orbit 
of the binary star / Moon. The latter interaction becomes important only when the 
protostellar disk fills the effective Roche radius of the star-disk system, and so may not 
be the dominant driver of orbital evolution. Orbital expansion may also depend on the 
eccentricity of this initial orbit. Note that binaries born from disk fragmentation may 
well be initially eccentric due to the m=l mode perturbations that can seed the disk 
dShu et al.|1990l|Krumholz et al.|2007HKratter et al. 2 010a) . 



3.2. Dynamical Interactions in Few Body Systems 

As mentioned above in the context of the competitive accretion models, dynamical 
interactions in small-n systems may be extremely important in reshaping the orbital 
distributions of young systems. One dominant dynamical mechanism for producing 
close binaries may be interactions among three or more stars. Binary-binary and single- 
binary encounters can lead to the ejection of (usually) the lowest mass star, leaving 
behind a tighter pair. Similarly, triple and higher multiplicity systems may form or 
evolve into unstable hierarchical configurations, where their relative orbits interact and 
trigger chaotic orbital evolution (Mardlin g & Aarseth|2001 1. Chaotic orbits can lead to 
ejections, tightening of the remaining binary, partner exchanges, eccentricity excitation, 
and in some cases, stellar collisions. These dynamical interactions are often seen in the 
competitive accretion simulations, and reproduce some features of the observed binary 
distribution of very low mass objects ( |Dupuy & Liu|20lT i. 

An even broader range of dynamical interactions may occur throughout the life- 
time of the binary, particular in dense stellar clusters. For a more detailed analysis of the 
different processes at work, and population synthesis calculations of interaction rates, 
see for example |Portegies Zwart et aL] ( |1997| ), [Pooley et al.| ( |2003| ), and [Fregeau et al] 
( [20041 ) 



We now explore the observational evidence for these theories in both the pre-main 
sequence and field populations. 



4. Linking Theory with Observations 

Advances in numerical techniques mean that each of the formation and evolutionary 
processes described above has been simulated at high resolution (see references above). 
Though numerical artifacts remain, many of the mechanisms cannot be ruled out on ba- 
sic physical grounds. Instead they must be evaluated based on their ability to reproduce 
the observables of the star formation process. 

An unequivocal determination of the binary formation process remains elusive; it 
would require that we know the conditions in a given cloud and core over its entire 
formation timescale. We have made great strides in characterizing the different phases 
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of core formation (Enoch et al. 2009 1, but we do not yet see binary formation in progress 



dSchnee et al.TIOTO 



We must, however, exercise some caution when assessing theories of binary for- 
mation that reproduce observed trends without accounting for the underlying physics. 
While the effort to create models that reproduce observations is important, we should 
not neglect detailed modeling of molecular clouds and prestellar cores from which bi- 
naries form. 



4.1. Observed Populations of Binaries and Multiples 

The field population of binary stars has a separation distribution that is log-normal 
for solar type stars, and bimodal and log-flat (/(log(P) oc const) for more massive 
systems (Ragh avan et al.|2010t|Sana & Evans |20l"0t [Mason et al.|2009] ). For the former 
population, the separation distribution is log normal, with a peak at fi = 5.03 days, and 
a standard deviation of criogp = 2.28. This leaves a sizable population within ~ 10 
AU that may be brought even closer through binary stellar evolution processes ( |Iben^ 
Tutukov 1985 1. The high mass stars include a large population of short period systems: 
as many as 60% of the systems have orbital periods less than 10 days ( [Sana & Evans 
2010|l. This overabundance of short period systems may be influenced by selection 



effects, since it is difficult to find lower mass companions to massive stars at distance 
of lO's-lOO's of AU dMason et al.|2009] ). 

While the mass ratio distribution of the low mass stars is gaussian (jj. ~ 0.6, cr ~ 
0.1), for more massive stars it is better represented by a power law f{q) oc q~^-'^, where 
q = m\lm2 ( Kobulnicky & Fryer 2007) . In both cases there is a bias towards more 
equal mass binaries at small separations. 

In contrast, recent work on the young cluster Taurus-Auriga paints a different pic- 
ture of the distribution ( Kraus et al.|[201 1) . These authors find a much higher binarity 
fraction of 73% even among low mass stars (< 2.5Mo), compared to 40-50% in the 
field. While the separation distribution is consistent with a log-normal, as in the field, 
they find that lower mass systems (M < 0.1 Mq) have a tighter separation distribution 
than their higher mass counterparts. Other young clusters show an enhanced binary 



fraction for low mass stars (Kraus et al. 2008 1, though Taurus is somewhat more ex- 
treme. 



4.2. Empirical Models for Turbulent Fragmentation 



Semi-empirical models such as Fisher (2004) have had some degree of success pro- 



ducing observed period and mass ratio distributions of low mass stars from the tur- 
bulent fragmentation mechanism. This work assumes that the components of binaries 
are drawn randomly from the initial mass function (IMF), and that there exists a con- 
stant efficiency factor with which mass and angular momentum are transferred from 
the core to the binary. They also draw eccentricities randomly from a thermal distri- 
bution, consistent with observations, in order to specify the orbit for a given binary 
angular momentum. A limitation of this work is that the physical process that gov- 
erns the distribution of mass and angular momentum between binary components is ad 
hoc, and tuned to reproduce the observations. While numerical work has shown that a 
multiple system can form out of a turbulent core, it remains to be shown in detail that 
the efficiency of angular momentum and mass removal scales as this model suggests. 
Moreover, while the final distribution of binary systems is consistent with the observa- 
tion distribution of low mass stars, it is not consistent with the orbital parameters and 
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mass ratios for substellar companions and more massive stars. Recent AO surveys by 



Metchev & Hillenbrand ( 2009 1 show that the distribution of substellar companions to 
solar-type stars is inconsistent with being randomly drawn from the IMF. There is also 
a dearth of substellar mass companions at large distances. On the opposite end of the 
mass spectrum there may also be inconsistencies. First, there is no explanation for the 



increase in binarity with core mass. Secondly, recent numerical work (Krumholz et al. 



2007 1 has shown that radiation from the first star can suppress turbulent fragmentation 
in massive cores because more massive protostars are luminous even at early times. 



4.3. Competitive Accretion and the IMF 

Similarly, the competitive accretion picture has had success explaining portions of the 
observed initial mass function and some binary statistics ( Moeckel & Bate|2010 l. How- 
ever, these models tend to overproduce the binaries containing the lowest mass stars and 



brown dwarfs, and have higher star formation efficiencies than observed (Price & Bate 



20091. Another criticism is that the relative velocities between accreting cores and 



background gas are so high that the accretion rate onto the cores should be far too low 
to produce more massive stars ( Krumholz et al.|[2006 ). In addition, the inter-core ve- 
locities appear to be much higher in the simulations than observed between actual star 
forming cores, calling into question whether or not there can be significant core-core 
interaction on the timescale of star formation (Kirk et al.|2010 Evans et al. 2009 ). 

The [Goodwin & Kroupa ( 2005[ ) picture does reproduce the observations for stars 
< IMo, but like the model of Fisher (12004) it has been tuned to do so. The mechanism 
which breaks some cores into two versus three stars, and the division of mass and 
angular momentum between them, is not specified in detail. 



4.4. Disk Fragmentation: a possible solution to the luminosity and angular mo- 
mentum problems 

Simulations now suggest that the disk scenario can produce a range of stellar sepa- 
rations and mass ratios, depending on primary mass, and the epoch of disk instabil- 
ity. Although these models do not make clear predictions for the final separation and 
mass ratio distribution of stars, the disk fragmentation scenario may be consistent with 



two observational puzzles in star formation: the so-called luminosity problem (Kenyon 
et al. I990j ) and the angular momentum problem (see Bodenheimer| 1 995" and references 
therein). 



4.4.1. The Luminosity Problem 

The luminosity problem refers to the discrepancy between the observed and predicted 
bolometric luminosity of young protostars based on their accretion luminosities. Given 
a star formation timescale and stellar mass, we can predict the time-averaged mass ac- 
cretion rate, and thus the corresponding luminosity. While uncertainties in the opacities 
and geometry of high extinction cores can complicate the interpretations, recent data 
from the Spitzer c2d survey confirms that most young sources are under-luminous com- 
pared to the expected luminosity from the time-averaged rate ( [Evans et^ 2009[ ). These 
authors suggest that variable, high accretion rates, particularly in the Class O/l phase, 
are all but required by the data. Because so few cores are observed to have luminosities 
consistent with the time averaged rate (a few 10~^Mo/yr) accretion must occur in short 
bursts (~ lO'^ years) well above this rate. However, a recent reanalysis has shown that 
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such extreme variability may not be required. Instead, OfFner & McKee (201 1 1 suggest 
that a time variable accretion rate can explain the bulk of the variability. 

For more massive protostars, high, variable accretion rates at early times are likely 
to produce fragmentation. However, such high mass cores are not typical among the 
sample described above, and thus the statistics on their observed accretion rates are 



limited (Klaassen et al.||20ll| ). Highly variable accretion from turbulent cores could 
further promote fragmentation if the disk is fed asymmetrically, which may also be 
consistent with observations (see e.g. Tobin et al. |2010[ and also Stamatellos et al. 



4.4.2. The Angular Momentum Problem 

The angular momentum problem refers to the fact that we do not understand how four 
orders of magnitude of angular momentum are repartitioned between the largest scales 
in Giant Molecular Clouds (GMCs) and the scales of binary orbits. By the time dense 
cores have formed, line-width size observations reveal that the typical specific angular 
momentum has dropped by two orders of magnitude ( Goodman et al.|1993 1. When we 
look at binary stars, whose formation is intimately tied to the core angular momentum, 
the average specific angular momentum can drop another two orders of magnitude. 

Magnetic fields may play an important role in removing angular momentum on 
the largest scales ( Li et al.|2004 i, although hydrodynamic turbulence can be significant 
as well ( Jappsen & Klessen||2004 1. Disk fragmentation may be important on smaller 
scales: because disk size is regulated by the gas reservoir's angular momentum, if the 
disk fragments to form a binary or multiple before all of the infalling gas and angular 
momentum has accreted onto the star-disk system, binaries can form with less orbital 
angular momentum than that contained in the core. 

4.4.3. Other Observational Support for Disk Fragmentation 

Several other trends seem roughly consistent with the disk scenario. Both disk frag- 
mentation and binary frequency increase with primary mass. At lower masses where 
disks are stabilized by stellar irradiation ( Matzner & Levin||2005 1, the turbulent cores 
out of which they form may be more susceptible to early fragmentation ([OfFner et al. 



20101 



The trend towards equal mass components at very high masses might be related to 



the fact that massive disks become unstable early on in their accretion histories (Pinson- 



neault & Stanek 2006). There are also hints in the high mass multiplicity statistics that 



mass ratio scales inversely with separation ( ,Mason et al.|2009| ). Although migration can 
alter the orbits, disk fragmentation suggests a possible correlation between secondary 
mass and separation: fragmentation typically occurs towards the outer edge of the disk 
which is larger at late times, when there is also less mass available for the secondary to 
accrete. 

Another piece of evidence favorable for the disk scenario is the observation that 
binaries separated by less than roughly 40 AU are more likely to have the stellar spins 
aligned with the binary orbit (Hale 1994[ ). Although this study does not account for 
higher order multiples, and has significant uncertainties in the measurements (e.g. sin 
i), this is suggestive that some objects were born in disks. Of course objects born from 
clouds with the same net rotation might also tend to be aligned. Recent observations of 
Herbig Ae/Be binaries also find alignment between the orbital plane of the binary and 
the accretion disks ( Wheelwright et al.|2011 1. 
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5. Conclusions 

Despite their ubiquity, we still lack a coherent theory of binary and multiple star forma- 
tion that can explain the enormous range of observed systems. This review provides a 
glimpse of the theoretical issues that must be resolved in order to make progress in the 
field. Interested readers should find ample references to explore these topics further. 

From a theoretical perspective, it remains difficult to predict the primordial dis- 
tributions in mass, separation, and eccentricity, of binaries and multiples. These dis- 
tributions must then be convolved with the known evolutionary processes in order to 
provide reliable initial conditions for binary stellar evolution calculations. At present, 
the observational data provides the best indication of the likely initial conditions for 
these models, despite lingering uncertainties regarding observational biases. 

Recent work shows that disk fragmentation requires careful consideration as one 
of the viable mechanisms, especially for stars more massive than our sun. In partic- 
ular, the combination of turbulent fragmentation across the mass spectrum with disk 
fragmentation added to the distribution at higher masses, might explain the increase in 
multiplicity fraction with stellar mass. 

Acknowledgments. The author would like to thank the conference organizers, and 
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